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In this work, we apply the vector meson dominance (VMD) model to extract the electromagnetic 
time-like form factor of the D* — > D transition combining the recent Belle data for e + e~ — > D* + D~ + 
c.c. and data for D* — > Z>y. Two solutions are obtained in the interpretation of the cross section 
lineshape: i) With a relatively large coupling for ipD* D determined by experiment, destructive 
interferences among those charmonium components are required to bring down the overall cross 
sections, and then account for the cross section lineshape. ii) With a relatively small value for the 
tpD*D coupling based on heavy quark theory, an apparent cross section deficit near threshold is 
observed, and contributions from other mechanisms are needed. It might imply the presence of an 
additional resonance X(3900). Meanwhile, we also point out that an enhancement like that could 
be produced by the D*D S + c.c. open channel effects. 

PACS numbers: 12.40.Vv, 13.40. Gp, 13.66.Bc 



I. INTRODUCTION 

The DD* + c.c. productions in e + e~ annihilation give access to the study of the time-like electromag- 
netic (EM) form factor of D* —¥ D^f* transition in the charmonium mass region. Their cross sections were 
measured recently by Belle [l[ and BABAR 0, and clear resonance structures were observed above the 
DD* or D*D threshold. In the real photon limit, the coupling form factor can be measured via D* — > Dj, 
which turns out to be an important decay mode for both the charged and neutral D* mesons [3j. In 
particular, it shows that the partial decay coupling for D*° — > D°j could be much larger than that for 
£)*± £)±^ This feature initiated a lot of efforts on understanding the D* — > D transition form factor. 

Our motivation in this work is to study the D* — > Dj* form factor in the time-like region with the 
help of the recent experimental data p], Q ■ We shall take into account the resonance contributions to 
the form factor by employing the extended vector meson dominance (VMD) model [1, [E] . To connect the 
real photon limit to the energy region above the D*D threshold, we also include the light vector meson 
contributions. In Ref. @, a VMD model was adopted for studying the D* — > £>y* form factor. However, 
due to lack of experimental information at that time, the authors assumed that the widths for all vector 
mesons apart from the ^(4040) (and beyond) are zero. This should be a too-rough approximation. As 
studied recently in Ref. [7[ , the width effects were found essentially important for understanding the cross 
section lineshape of e + e~ — > DD. 

Another useful and correlated channel is D* — > De + e~ , which probes the time-like form factor in small 
momentum squared region. However, due to the significant suppression of the EM vertex, branching ratio 
of this channel is expected to be very small and hard to measure. This branching ratio can be calculated 
in our model and serves as a prediction from theory. 

There is a great advantage for extracting the D* — > Dj* form factor in e + e~ — > D*D + c.c. Namely, 
there is only one Lorentz structure for the VVP coupling, where V and P stand for vector and pseu- 
doscalar meson fields, respectively. Therefore, all information about the transition mechanisms would 



* E-mail: yjzhang@ihep.ac.cn 
t E-mail: zhaoq@ihep.ac.cn 



2 



be contained in a single coupling form factor, which is a complex function of the photon's four-vector 
momentum squared. Our calculations will be compared with the Belle data for e + e~ — » D + D*~ +c.c. 

As follows, we first present the details of the VMD model in Sec. HU The numerical results will be 
given in Sec. HID Section [IV] is devoted to a summary and discussion. 



II. THE MODEL 



The typical effective Lagrangian for the j*D*D and j*D*D coupling can be written as: 

C = -ieg rD , B £ a0fiV d a A d"D* v D + h.c, (1) 

where A 13 is the vector meson and electromagnetic held, Sap^u is the antisymmetric tensor. With Eq. (JXJ) , 
the matrix element of e + e~ — > D*D in the one-photon approximation can be written as: 

T = e 2 v(k 2 )j a u(k 1 )^g 1 , D , B (s)e a f St _ lu p B p D ,e u , (2) 

where u(fci) and v{k2) are the Dirac spinors of the electron and positron, respectively; e" represents 
the Z)*-meson polarization vector, and g 1 * D * B {s) is the effective coupling form factor for the D* —> D 
transition. Note that the electron charge e has been isolated out in this definition. In the above equation, 
s = (ki + &2) 2 is the overall center mass energy, while p B and po* are the four-vector momenta of the 
final state D and D* meson. 





FIG. 1: Schematic diagrams based on the VMD model for e + e~ — > D* D + c.c. Diagram (a) is for the single 
photon approximation with an effective coupling P 7 *_d*_d(s), while (b) represents that the electromagnetic field is 
decomposed into the sum of vector meson fields. 



As shown in Fig. [Q with the VMD model H U we can decompose the electromagnetic current into 
a sum of all vector meson fields including both isospin-0 and isospin-I components. The V^* effective 
coupling can be written as: 

Cvy = J2 e -^V^, (3) 
where V M (= p, ui, <fi, J/i/j...) is the vector meson field, and eMy/fy is the photon-vector-meson coupling 
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TABLE I: Resonance parameters of the vector mesons adopted in this study. They are taken from PDG @]. 





p(770) 


/o(1450) 


p(1700) 


u(782) 


cj(1420) 


cj(1680) 


J/4* 


V(3686) 


V>(3770) 


1^(4040) 


M v (GeV) 


0.774 


1.465 


1.720 


0.783 


1.45 


1.62 


3.097 


3.686 


3.773 


4.039 


IV(MeV) 


149.4 


400.0 


250.0 


8.5 


200.0 


250.0 


9.32 xl0 _ii 


0.317 


27.3 


80 


T ee (keV) 


7.04 






0.6 


0.46 


0.8 


5.55 


2.38 


0.265 


0.86 



constant. Setting m e 



0, e//v can he extracted from the partial decay width IV_> e + e - 

i 

3rv_s. e + e - 



by: 



e 



2a e \p e 



(4) 



where \pl\ is the electron three- vector momentum in the vector meson rest frame, and a e is the fine- 
structure constant. 

The following effective Lagrangians are required for vector meson couplings to the D meson pair and 
D*D: 



£-VDD 
£-VD*D 



gvDSiDdnD-dnDDyV, 
-ig VD ,D^ v d a V^d^D* v D + h.c. 



The effective coupling <? 7 *£)»£i(s) can then be expressed in a general form: 



Ml 



'9VD'D' 



where IV is the total decay width of the vector meson. The total cross section for e + e 
thus reads 



D*+D~ 



(5) 



(0) 



- c.c. 



a(e + t 



-> D* + D~ +c.c. 



3 s 3 / 2 ' 
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fv S- 



My 
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whereg VD , £, = g VD ,+ D - — gvD*-D+- I n this paper, our definition for gvn±£> 1S different from that in 
Ref. Q- Namely, it does not include the charge conjugate coupling. In Ref. |3, gvD*D = V%gvD'+D- = 
V2gvD*-D+- Thus, Eq. Q has a factor of 2 different from Eq. (12) in Ref. 0. 

In the following calculation, we mainly consider the contributions from p, w, J/tp, and their radial 
excitation states, the contributions from the 0-mesons are dropped because the g^u-Q couplings are 
strongly suppressed by the Okubo-Zweig-Iizuka (OZI) rule. The g-fD*D couplings are also suppressed by 
the OZI rule. Moreover, the T states are far away from the D*D threshold. Thus, their contributions 
can be safely neglected. Parameters for the vector mesons are listed in Table |U 

The asymptotic behavior of g 1 , D ,jj{s) has been discussed in Ref. @. As s — > oo, form factor <7 7 »£)»£>(s) 
must decrease at least as s -2 to avoid the violation of unitary. As a consequence, the following relations 
are obtained for the isospin-1 and the isospin-0 components: 



Pi 



— f — 9 Pi D*D — ) 

Jpi 



and 



M 2 

IVI V(I=0) 



V(I=0) J 



V(I=0) 



gv(i=o)D*D — o 



(8) 



(9) 



As discussed in Ref. [f| , the above asymptotic relation implies that at least two p mesons are needed in 
the VMD model. 



4 



We also adopt the following relations given by the SU(3) quark model: 

fu t - 3/ w , m 2 p . ~ ml. , (10) 

where the factor 3 can be tested well by the partial decay widths for uti and pi via Eq. ((4]). With the 
flavor symmetry, we also have: 

SuiiD' + D- = -9 Pi D* + D-- (11) 

It is worth noting that the above relations, i.e. Eqs. ([8]), (fT0|) and (fTTj). imply a negligible contribution 
from the p and u) mesons in the production of D*D + c.c. pairs in e + e™ annihilation, although they are 
dominant in the D* radiative decays. 

The strong coupling gvD'D f° r those charmonium states below the D*D threshold cannot be directly 
extracted from the experimental data, such as gj/^o-D- Their coupling values generally have large 
discrepancies in different models. The relation between gvD*D an d 9vdd can b e parameterized by: 

dVD'D = 9vdd x 9v, (12) 

where gy is a parameter with inverse of mass dimension. In the heavy quark mass limit [8[, one has 
gv = \j Jmrywin ~ 0.52 GeV -1 . In contrast, the relativistic potential model of Ref. [9| gives: 

where cq is the heavy quark charge and e q is the light quark charges, and the expressions of Aq and 
A 9 can be found in Ref. @. It should be mentioned that as pointed out in Ref. [IoT ]. Equation (TIB"]) is 
a general consequence of decomposing the electromagnetic current into a heavy and light part in the 
radiative decay D* — > Z>y. In Ref. [Tl(, the value gj/^ ~ 1.0 GeV -1 is extracted, and QCD sum rules 
give gj/j, = 0.69 ± 0.14 GeV -1 pi]. 

In the numerical calculation, we neglect the coupling differences between the charge and neutral chan- 
nels for gvD'D an d 9vdd- The coupling constants g^ijjjj and <?,/, (3770)0.0 have been discussed in our 
previous work Q, and 3^(3770)00 — 12-7 is extracted from the experimental result Q by the effective 
Lagrangian approach. The coupling g^uT) — 9-05 is determined by fitting the lineshape of e + e~ — > DD 
process. We adopt gj/^oD = 7-44 from Ref. which is obtained by the VMD model. 

We must note that the cross section measurement gives access to the absolute value of |3 7 *_D*r>( s )l- But 
9-j'D*d( s ) i s a complex function of s in the time-like region. Moreover, the prescription at the hadronic 
level will introduce a phase factor e 1 ^ to each resonance amplitude. These phase angles, apart from an 
overall phase, can be determined by fitting the Belle data for e + e~ — » D* + D~ + c.c. [l(. 

As one can see that the Belle data cover a rather high ^/s region, it is natural to anticipate that the 
low t/s form factor would be less sensitive to the data constraints. Taking into account this, we include 
the real photon data for D* — > Dj in the numerical fitting. In the real photon limit, it is rather direct 
to obtain the partial width of D* — > D7 by 

r(zr D 7 ) = ?fg 2 lD * D (0)q 3 , (14) 

where q is the photon energy in the D* rest frame. Again, <7 7 o*o(0) can be expressed as Eq. (J6j) . This 
would provide constraints on the parameters from light vector meson components. 

The decay of D* —> De + e~ also gives access to the transition form factor at low y/s. The matrix 
element of the D* — > De + e~ decay via single photon transition can be written as 

e 2 

T = 9-y* d* D( k2 )-^ £ af!^PDPD' €L 'u(k 1 )j a v(k 2 ), (15) 

where k 2 is the invariant mass of the lepton pair. The formula of the differential probability can be 
described by the following expression: 



dT = a 2 e \ gj , D , D ( k 2 )\ i; , w o , 
dk 2 dQ 2 487rm D 



48 [k +2Q 2m D 2m D , 
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8m 2 e Q 2 



(16) 



where Q 2 is defined as Q 2 — (k 2 +Pd) 2 (or Q 2 = (k% +pn) 2 )- 
suppressed by an additional EM coupling in respect of D* - 
measure this branching ratio in experiment. 



We should note that this process is strongly 
-> Dj. Therefore, it is relatively difficult to 



III. NUMERICAL RESULTS 



Now, we switch to the details of the numerical fitting. First, we give a brief discussion about the fitting 
scheme: 

(i) Since only the relative phase can be measured in experiment, we set (jjj/^ — 0°, and then the other 
phase angles are defined in respect of <f>j/^. Meanwhile, since the cross sections are not sensitive to the 
light vector meson contributions, the relative phases <p Pi and tf> Ui are set the same and denoted by (fi^y 
in the fitting. 

(ii) As shown by the cross sections around 4.2 GeV, there is no clear evidence for an enhancement due 
to the presence of a resonance. Therefore, the data for e + e~ D*D + c.c. cannot constrain 0^(4260) a t 
all. For simplicity, we set 0^(4260) — 0°. 

(iii) For -0(4040) and F(4260), gvD'D are n °t clear. Especially, the y^- for V(4260) is also unavailable. 
In the numerical fitting, these two couplings are always combined. Thus, we define g^ = y^- x gyD'D 
for -0(4040) and Y(4260). In total, the fitting parameters include the relative phases 4>lv, 4>4>', 0v>(377o); 

</y>(4040) , and the couplings ^f 404 o) and ^(Leo)- 

As discussed previously, the charmed meson couplings to the light mesons are obtained in the chiral 
and heavy quark limits [14]: g PlD *D = s/2\m Pi /f T , with / w = 132 MeV, and A = 0.56 GeV -1 These 
couplings contain uncertainties arising from gy as illustrated by Eq. (|13|) . Also, since the constraints on 
the light vector mesons are rather weak in the data for e + e~ — > D*D + c.c, we thus include the data for 
D* — > Dj to constrain couplings g p and g u . 

For charmonium coupling to the charmed mesons, the following relation is assumed: 

SV> = 9J/ip — 9if>' — Si/>(3770) — ffi>(4040) ) (17) 

which can be determined by 

r(0(4O4O) -> DD) |pl| 3 



r(^(4040) -> D*D + c.c.) 4M 2 (4040) |p 2 | 



(18) 



where p[ and p2 are the three momenta of the final charmed mesons in 0(4040) — > DD and -0(4040) — > 
D*D + c.c, respectively. Hence, given the experimental data [||, 

r(0(4O4o M W) =Q05±0 03; (19) 



r(^(4040) ->• D* D° + c.c.) 



we have = 1.73 ± 0.52 GeV 1 by taking the average value corresponding to the datum bound. This 
value appears to be larger than gy ~ 0.52 GeV -1 extracted by Ref. || and gy ~ 1.0 GeV -1 by 
Ref. [Ill ] . In order to examine the impact of the uncertainties due to gy, we shall fix g^ = 1.73 GeV -1 
and 0.52 GeV -1 , respectively, in the numerical fitting. 
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A. With relatively large g^, determined by experimental data 



In Table HH all the fitted parameters are listed. It shows that the coupling g^ i2 &o) nas a l ar S e 
uncertainty, which reflects the negligible role played by V(4260) in the fitting. Further experimental data 
with high accuracy are needed to extract its resonance parameters. 



TABLE II: Model parameters obtained from the x 2 minimization fitting with = 1.73 GeV 1 . Coupling gy 
is defined by g^ = x gvD-D- The phase angles are in radian. The reduced x 2 is x 2 /d.o.i = 41.2/51. 



parameter 


4*LV 


<t>i>> 


0V(377O) 


<0/>(4O4O) 


fli/'(4040) 


5y(4260) 




— 0.41 ± 0.17 


2.36 ±0.32 


4.54 ± 0.22 


-0.68 ±0.21 


0.37 ±0.05 


0.02 ±0.03 



With the help of Eq. (g]) and the fitted value for g^( 4040 ), we obtain 3^,(4o4o)_d*d — 0-74 ±0.1 GeV . 
Consequently, the partial decay width of -0(4040) — > D*D ± c.c. can be accessed: 

r(>(4040) ->■ D* + D~ +c.c.) = 5.1 ± 1.0 MeV, 

r(i/i(4040) -> D*°D° ± c.c.) = 5.5 ± 1.1 MeV. (20) 

An interesting result from this fitting is that, although gyo'D bares large uncertainties, the excita- 
tions of the charmonium states J/0, 0' '> & n d "0(3770), and their interferences play a major role on the 
interpretation of the effective coupling g 1 *£,»£ ) or g^j^D 0- This feature can be seen more clearly via 
the fitted cross sections. 

In Fig. HJa), the total cross section and cross sections for exclusive resonances are plotted. We do not 
show the curve of V(4260) since its contribution is negligibly small. Interestingly, other charmonia, such 
as J/tp, 0', and 0(3770), have large exclusive cross sections. In particular, the cross section for the 0' 
excitation over-shoots the data apparently, and cancellations among these three amplitudes are required 
to reproduce the lineshape of the cross sections. 




FIG. 2: The Belle data for the e + e - — > D* + D~ ± c.c. cross section [lj are fitted by the x? minimization method 
with different values for g^. Panel (a) is obtained with g^ — 1.73 GeV -1 , and panel (b) with g^ — 0.52 GeV -1 . 
The dotted, dashed, short-dashed and short-dotted lines are for exclusive contributions from the light-vector 
mesons (LV), J/ijj, ip' and 0(3770), respectively. In panel (a), the solid line represents the overall results, while 
the dash-dotted line if for exclusive contribution from 0(4040). In panel (b), the solid, dash-dot-dotted lines are 
for the overall results from two fitting schemes, i.e. Scheme-A and Scheme-B, respectively, while the dash-dotted 
and short-dash-dotted are for contributions from t/>(4040) in these two schemes. 
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We can then extract the form-factor \g 7 *D*+D- I m the whole time-like region. In Fig. EJa), we first 
look at the region around the D*D threshold up to 5.0 GeV. As demonstrated by the solid line, the data 
can be described perfectly. We also include those two empirical fits presented in Ref. Q as a comparison. 
The dashed line is generated by fitting the form factor data with an exponential function, i.e. with form 
factor one (FF-I): 



V2g^, D ,+ D -(s) = ffiexp [-(s - (m D + m D *) 2 )/t 1 ] + g , 



(21) 



where x — corresponds to the D*D + c.c. threshold, and g\ 1 t\, and go are fitting parameters. The 
dotted line is given by fitting the data with a single resonance: 



V2g 



7 *D* + D 



■(«) 



s — m x + imx^x 



(22) 



with a background term b\. The parameter b can be regarded as the product of the 7* A coupling 
and XDD* coupling. This parametrization agrees with the data at higher energies, but drops at the 
threshold. In the above two equations, a factor \[2 has been included for the change of conventions here. 
All the parameters have been given in Ref. 0. 

In Fig. |U[a), we plot the v^-dependence of the form factor g 1 -£,-o{s) for both charged and neutral 
channel in the whole time-like region up to 5.0 GeV. Since the coupling g^o'D (and gj/^o'D) nas t ne 
same value in these two channels, these two form factors converge to each other around the threshold 
region. In the low-i/s region, the discrepancy arises from the total width difference between D* ± — > D^j 
and D*° — > D 0r y. For the latter, the experimental data only give an upper limit, i.e. T(D*° — > D°j) < 2.1 
MeV 0]. This corresponds to |5 7 _do»d° (0)| < 11-3 GeV -1 , which has not been marked in the figure. 

In Fig- Hfb), the form factor in the small ^/s region is plotted in association with its real and imaginary 
part. The resonance structures from p(770), p(1450) and p(1700) are distinguishable. However, it should 
be cautioned that this kinematic region would suffer from a lack of information about the light vector 
meson couplings to D*D. Our model bridges the real photon form factor with the high-y / s one, but 
inevitably leaves the middle range with large uncertainties. We drop the high-y / s part between 2.5 GeV 
and 4.5 GeV since the structure of the real and imaginary part appears trivially as either very narrow 
peaks or very narrow dips. 




FIG. 3: The y's-dependence of form factor |.9 7 *_d* + _d- ( s )I extracted from the fitting results. In panel (a), the 
solid line is obtained from the fitting results with = 1.73 GeV -1 , while the dashed and dotted lines are for 
FF-I and FF-II results from Ref. [7(. In panel(b), the dashed and the dotted lines are obtained from the fitting 
results with — 0.52 GeV -1 for Scheme-A and B, respectively. 

The following points are advocated to understanding this fitting results: 
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(i) Numerically, the large contributions from ip' are due to a relatively large g^ as shown by Eq. (fl2")) . 
In Ref. 0, we showed that the coupling j^^g can be well-constrained by the cross section lineshape of 
e + e - — > DD. Thus, the coupling g^ijyjj is actually enhanced by g^ via Eq. (fT2")l . 

(ii) The relative phases appear to be sensitive to the cross section lineshape of e + e~ —> D*D, and can- 
cellations among the dominant amplitudes seem to be inevitable. Such interferences generally would affect 
the extraction of resonance parameters. Because of this, it is desirable to have a precise measurement of 
the lineshape of e + e~ — > D*D + c.c 
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FIG. 4: The y^s-dependence of form factor |s 7 »_d»_d(s)| extracted from the fitting results with — 1.73 GeV -1 . In 
panel (a), the dashed and dotted line stand for the form factors for the charged and neutral channel, respectively. 
In panel (b), the solid line is the form factor \g, f , D *+ D - (s)| for the charged channel, while the dashed and dotted 
line denote respectively the real and imaginary part of g 7 *£>* + zj- ( s )- 

In order to clarify the interferences among the charmonium states, we identify the exclusive contribu- 
tions from light vector mesons and charmonium states in D*^ — > I? ± 7, where we expect that the light 
vector mesons should play a significant role. We first inspect the separated contributions to the partial 
width of D*^ _D ± 7 from the light vector mesons and charmonia, and the results are as follows: 

rn/(L>* ± -> D ± 1 ) = 2.21 keV, 
r^,(Z>* ± -> D ± 1 ) = 1.61 keV, (23) 

where T^y and are decay widths contributed by the light vector mesons, i.e. p, to etc, and the 
charmonia, i.e. J/ip, 4>' etc- Interestingly, it shows that, although the light vector mesons play a dominant 
role, contributions from the charmonium states are still sizeable. In comparison with the experimental 
result r(L>* ± D± 7 ) = 1.54 keV @, one can see that a destructive interference between these two 
components is required. 

Such destructive phases are also present within the charmonium states. As follows, we list the exclusive 
contributions from J/ip, ip', and ■)/'(3770): 

Tj/^D^-tB^j) = 8.07 keV, 
iyCD^-^D^) = 5.11 keV, 
IVpttojP** ^D ± 1 ) = 0.93 keV. (24) 



Other charmonium exclusive contributions are negligibly small. So we do not list them here. 
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B. With smaller g^, by heavy quark theory 

As pointed out earlier, the relatively large contributions from the charmonium states, in particular, 
J/ip and jp' , are mainly due to the larger value of cu in comparison with those given by the heavy 
quark theory Q and the relativistic potential model [a] . In order to examine the impact from a possible 
overestimate of the coupling <jy, due to the experimental uncertainties @, we perform another fit of the 
e + e~ — > D* + D~ + c.c. cross sections adopting g^ = 0.52 GeV -1 §]. 

Two fitting schemes are considered. In Scheme- A, we fix g^ = 0.52 GeV -1 and fit the data with the 
same parameters. The fitted parameters are listed in Table IIVI Relatively large x 2 is found and the 
fit cannot account for the cross section lineshape near threshold as shown by the solid line in Fig. IS^b). 
The major deviations occur near threshold where a cross section deficit is revealed. The relatively small 
g^, also leads to suppressed contributions from J/ip, ip' and ^(3770). In contrast with the fitting of 
Subsection IIII Al this scheme requires constructive phases among the charmonium amplitudes. 

In Schcmc-B, we still fix g^ — 0.52 GeV" 1 , and introduce contributions from an additional resonance 
AT (3900) in order to overcome the cross section deficit near threshold. When we introduce the X(3900) 
resonance, it is difficult to obtain the fitting solution if we leave the resonance parameters free, i.e. the 
mass and total width. We then take mx — 3.9 GeV and Yx — 89.8 MeV, which were extracted from 
e + e~ — > DD in Ref. Q. as an input. The other fitted parameters are listed in Table ITVl and relatively 
smaller \ 2 1S found. In Fig. [U[b), the results for Scheme- A and B are compared with each other. It 
shows that contributions from the -0(4040) are relatively suppressed due to the presence of the V(3900). 
Some interfering effects for the V(4260) are also observed in Scheme-A, but not in Scheme-B. As shown 
in Table HVl the combined coupling for the V(4260) still has large uncertainties. 



TABLE III: Model parameters obtained from the \ 2 minimization fitting with = 0.52 GeV 1 . We fix 
mx = 3.9 GeV and Tx = 89.8 MeV, which are from Ref. 0. The phase angles are in radian. 



parameter 


4>LV 


</V 


0^(3770) 


^(4040) 


ffj/;(4040) 


yy(4260) 


ax(3900) 


^X (3900) 


X 2 /d.o.f 


Scheme A 


5.72 ±0.63 


3.56 ±0.10 


3.69 ±0.22 


0.0 ±0.11 


0.58 ± 0.05 


0.08 ±0.03 






63.9/51 


Scheme B 


5.69 ±0.66 


4.21 ±0.43 


2.23 ±0.35 


5.62 ±0.27 


0.35 ±0.07 


0.01 ±0.04 


1.89 ± 0.44 


3.03 ±0.26 


38.8/49 



In Fig.^b), the form factors extracted from Scheme-A and B are compared with the data. Again, we 
see that a smaller value for g^, cannot account for the form factor near threshold. Although the inclusion 
of an additional resonance A" (3900) can optimize the description, the form factor appears to drop quickly 
in the subthreshold region. This will lead to different trends of the form factor in the middle range of 
kinematics, i.e. yfs is between 2 ~ 3.9 GeV. 

The results with a smaller value for gy, seem to be different from those with a larger value. The 
following points can be learned and conjectured: 

i) The resonance parameters extracted from the cross sections would contain uncertainties inevitably 
due to the uncertainty with g^. 

ii) The need of AT (3900) with a small <jy, would bring questions on the underlying physics. On the 
one hand, if such a resonance indeed exists, problem will arise from how to organize it within the quark 
model framework. The systematic study of quark potential model seems not to have a place for this state 
below 4.2 GeV PJ]. 

iii) To void the conflicts with the quark model, one possibility for such a structure would be due to the 
open channel effects of D*D S + c.c. Its threshold is 4.08 GeV, which is not far away from the D*D + c.c. 
threshold. The final state interaction D*D — > DD* via kaon exchange can produce a resonance- like 
enhancement near threshold. This is similar to the mechanism discussed in Ref. [jj. Also, justification 
of such a possibility would need the data for e + e~ — > D*D S + c.c, which unfortunately are unavailable. 

iv) As a comparison with the results from Subsection IIII Al we also list the exclusive contributions of 
the charmonium components to the radiative width in Table IIVI Again, we see that with the smaller g^ , 
the charmonium contributions are also strongly suppressed in the real photon form factor. Nevertheless, 
a destructive phase exists between the J/ip and ip' amplitudes. Note that the contributions from the 
light vector mesons are unchanged. 
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TABLE IV: Exclusive contributions to the D* ± — > D ±r y partial width from individual charmonium states with 
g^p = 0.52 GeV -1 . is the width given by a coherent sum of all the charmonium amplitudes. 



Partial width (keV) 






1V(3770) 


r^,(4040) 


r<0 


Scheme-A 


0.72 


0.46 


0.08 


0.01 


0.19 


Scheme-B 


0.72 


0.46 


0.08 


0.00 


0.31 



C. Predictions for D* ~ > De + e 



The ambiguity with the charmonium couplings will not affect the calculations for D* — > De + e - as 
long as the radiative decay D* — > D"f is fixed. This is simply because the mass of the electron (positron) 
is very small, and y/s of the virtual photon is close to the real photon limit. 

With the form factor determined in Subsection IIIIA1 we can predict the partial decay width of D* — > 
De + r with the help of Eq. ([T6|): 

rCD** D ± e + e-) = 9.95 eV. (25) 
It corresponds to a branching ratio, i?i?(Z?* ± — > D ± e + e~) = 1.04 x 10 -4 , and should be accessible in 



experiment, e.g. at BES-III 17]. This decay channel may provide some further constraints on the role 



played by the light vector mesons and charmonium states. 



IV. SUMMARY AND DISCUSSION 



In this work, we study the D* D + c.c. production in e + e - annihilation from the threshold to y/s ~ 5.0 
GeV in the VMD model. The recent experimental data from Belle [l[ allow us to extract the form factor 
9~{*d»d{ s ) m the time-like region. 

Due to the uncertainties with the charmonium couplings to D* D + c.c, we find two different solutions 
in the interpretation of the experimental data: i) With a relatively large coupling for i/jD*D, significantly 
large contributions from individual charmonium states are found. Destructive interferences among those 
charmonium components are hence required to bring down the overall cross sections, and then account 
for the cross section lineshape. ii) With a relatively small value for the ipD*D coupling, i.e. = 
0.52 GeV -1 Q, we find an apparent cross section deficit near threshold, and contributions from other 
mechanisms are needed. 

We also try to fit the e + e - — > D* + D~ + c.c. with = 0.69 GeV -1 , which is obtained by QCD sum 
rules [l2| ■ It shows that a cross section deficit still exists near threshold, although is becomes smaller 
than that with g^ = 0.52 GeV -1 , and x 2 /d.o.f = 54.6/51 is found. This suggests that a better fit of 
the near-threshold cross section would require a relatively large value for g^. The consequence, however, 
is that destructive inferences among the resonances beyond the threshold region would be expected. In 
case that g^ has a relatively small value, such a cross section deficit, on the one hand, might imply the 
presence of an additional resonance X(3900). On the other hand, we point out that an enhancement like 
that could be produced by the D*D S + c.c. open channel effects, and further experimental data will be 
able to clarify this issue. 

InRef. 0, 9sr* d*+ d~ ( s ) is extracted by two simple functions from which the behavior of the form factor 
below the D* D threshold is, in principle, unknown. In this work, our model provides a description of the 
form factor 3 7 *d»£)(s) in the interplay region between the real photon energy and D*D threshold, and 
some insights into the evolution of the vector meson contributions can be gained. Since this is a kinematic 
region which cannot be directly accessed by experiment, we expect further investigation of alternative 
processes would provide a test of our model, and more information on the underlying dynamics could 
be extracted. In particular, the partial decay width for D* — > De + e~ is predicted in this framework. 
With a branching ratio of Si?(£>* ± — > D ± e + e - ) = 1.04 x 10 -4 , this channel can be measured by BES- 
III in experiment. Our results also support such an ideal that the enhancement around 3.9 GeV in 
e + e - -t DD [H] is caused by the open D*D threshold effects 0- 
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It is worth mentioning a recent study of multiple solutions in extracting physics information from 
experiment (l9| . Since it is not possible that all the numerical solutions are correct, a conjecture of 
selecting the physical solution is that the physical solution would correspond to the minimal magnitudes 
of the amplitudes. Further theoretical investigation of the ipD*D couplings and experimental information 
for e + e _ — > D*D s +c.c. may help disentangle the underlying mechanism and also test the idea of Ref. (l9j . 
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